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Abstract 

The first-order reversal curve (FORC) method for analysis of sys- 
tems undergoing hysteresis is applied to dynamical models of electro- 
chemical adsorption. In this setting, the method can not only differ- 
entiate between discontinuous and continuous phase transitions, but 
can also quite accurately recover equilibrium behavior from dynamic 
analysis for systems with a continuous phase transition. Discontinu- 
ous and continuous phase transitions in a two-dimensional lattice-gas 
model are compared using the FORC method. The FORC diagram for 
a discontinuous phase transition is characterized by a negative (un- 
stable) region separating two positive (stable) regions, while such a 
negative region does not exist for continuous phase transitions. Ex- 
perimental data for FORC analysis could easily be obtained by simple 
reprogramming of a potentiostat designed for cyclic-voltammetry ex- 
periments. 

Keywords: First-order Reversal Curve; Hysteresis; Continuous phase 
transition; Discontinuous phase transition; Lattice-gas model; Monte Carlo 
simulation; Cyclic-voltammetry experiments. 
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1 Introduction 



Recent technological developments in electrochemical deposition have made 
possible experimental studies of atomic-scale dynamics jT]. It is therefore 
now both timely and important to develop new computational methods for 
the analysis of experimental adsorption dynamics. In this paper we apply 
one such analysis technique, the first-order reversal curve (FORC) method, to 
analyze model systems with continuous and discontinuous phase transitions. 
We propose that the method can be a useful new experimental tool in surface 
electro chemistry. 

The FORC method was originally conceived in connection with the 
Preisach model of magnetic hysteresis. It has since been applied to a vari- 
ety of magnetic systems, ranging from magnetic recording media and nanos- 
tructures to geomagnetic compounds, undergoing rate-independent (i.e., very 
slow) magnetization reversal [3]. Recently, there have also been several FORC 
studies of rate- dependent reversal jH El El- Here we introduce and apply 
the FORC method in an electrochemical context. For completeness, a brief 
translation to magnetic language is found in the Appendix. 

We apply FORC analysis to rate-dependent adsorption in two-dimensional 
lattice-gas models of electrochemical deposition. Specifically, we study a 
lattice-gas model with attractive nearest-neighbor interactions (a simple model 
of underpotential deposition, UPD), being driven across its discontinuous 
phase transition by a time- varying electrochemical potential. In addition, we 
consider a lattice-gas model with repulsive lateral interactions and nearest- 
neighbor exclusion (similar to the model of halide adsorption on Ag(lOO), 
described in Refs. 13 |Hl El CDl ) , being similarly driven across its continuous 
phase transition. 

The rest of this paper is organized as follows. In Sec.|21the FORC method 
is explained. The model used for both systems with continuous and discon- 
tinuous transitions is briefly discussed in Sec. El In Sec. jH the dynamics 
of systems with a discontinuous phase transition are studied using Kinetic 
Monte Carlo (KMC) simulations, as well as a mean-field model. The dy- 
namics of systems with a continuous phase transition are studied in Sec. El 
Finally, a comparison between the two kinds of phase transitions and our 
conclusions are presented in Sec. El 
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Figure 1: Schematic diagram of a FORC curve. 



2 The FORC Method 

For an electrochemical adsorption system, the FORC method consists of sat- 
urating the adsorbate coverage 6' in a strong positive (for anions; negative for 
cations) electrochemical potential (proportional to the electrode potential) 
and, in each case starting from saturation, decreasing the potential to a series 
of progressively more negative "reversal potentials" fir (Fig-IH). Subsequently, 
the potential is increased back to the saturating value PJ. It is thus a simple 
generalization of the standard cyclic voltammetry (CV) method, in which 
the negative return potential is decreased for each cycle. This produces a 
family of FORCs, 6{jlr,P'), where 6 is the adsorbate coverage, and where jl 
is the instantaneous potential during the increase back toward saturation. 
Although we shall not discuss this further here, it is of course also possible 
to fix the negative limiting electrode potential and change the positive return 
potential from cycle to cycle. 

It is further useful to calculate the FORC distribution. 



2 d^r <9/i ' 

which measures the sensitivity of the dynamics to the progress of reversal 
along the major loop.^ The FORC distribution is usually displayed as a 
contour plot called a 'FORC diagram.' A positive value of p indicates that 

^Note that to normalize the FORC distribution, the term ^5{fl ~ Mr ) | 
must be added to Eq. JI]. Here we consider the distribution only away from the line 
p, = jlr- The additional term could be found from the major loop. 
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the corresponding reversal curves are converging with increasing /i, while a 
negative value indicates divergence. Some preliminary results of this work 
have been submitted for publication elsewhere fT^ . 



3 Model 

KMC simulations of lattice-gas models, where a Monte Carlo (MC) step 
corresponds to an attempt to cross a free-energy barrier, have been used to 
simulate the kinetics of electrochemical systems with discontinuous El 
IT^ ITlj or continuous [71 phase transitions in two dimensions. The energy 
associated with a lattice-gas configuration is described by the grand-canonical 
effective Hamiltonian for an L x L square system of adsorption sites, 

= - ^ 4>ijCiCj - fi'^Ci , (2) 

i<j i=l 

where J2i<j is a sum over all pairs of sites, (pij are the lateral interaction 
energies between particles on the ith and jth sites measured in meV/pair, 
and yU is the electrochemical potential measured in meV/atom. The local 
occupation variables Cj can take the values 1 or 0, depending on whether site 
i is occupied by an ion (1) or solvated (0). The sign convention is chosen such 
that p, > favors adsorption, and negative values of 0jj denote repulsion 
while positive values denote attraction between adsorbate particles on the 
surface. In addition to adsorption/desorption steps, we include diffusion 
steps with a comparable free-energy barrier [7]. 

In each time step of the KMC simulation, an adsorption site is chosen 
at random and the transition rates from the present configuration to a set 
of new configurations (desorption, diffusion) are calculated. A weighted list 
for accepting each of these moves is constructed using Eq. ^ below, to 
calculate the probabilities i?(F|I) of the individual moves between the initial 
state I and final state F. The probability for the system to stay in the initial 
configuration is consequently -R(I|I) = 1— Sf^i-R(F|I) [71IH]. 

Using a thermally activated, stochastic barrier-hopping picture, the en- 
ergy of the transition state for a microscopic change from an initial state 
I to a final state F is approximated by the symmetric Butler- Volmer for- 
mula [iniiiiiiiHi 

Ct.=^ + A.. (3) 
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Here Ui and U-p are the energies of the initial and final states, respectively, 
T\ is the transition state for process A, and is a "bare" barrier associ- 
ated with process A. This process can here be either nearest-neighbor diffu- 
sion (Ann), next-nearest-neighbor diffusion (Annn), or adsorption/desorption 
(Aa/d)- The probability for a particle to make a transition from state I to 
state F is approximated by the one-step Arrhenius rate |TB1 IT7| ITH] 

-e.p ^ .e.p {-^) e.p {-"-0^) , (4) 

where u is the attempt frequency, which sets the overall timescale for the 
simulation. The electrochemical potential /i, which is proportional to the 
electrode potential, is increased monotonically, preventing the system from 
reaching equilibrium at the instantaneous value of fl. 

Independent of the diffusional degree of freedom, attractive interactions 
> 0) produce a discontinuous phase transition between a low-coverage 
phase at low /i, and a high-coverage phase at high fl. In contrast, repulsive 
interactions < 0) produce a continuous phase transition between a low- 
coverage disordered phase for low /i, and a high-coverage, ordered phase for 
high fi. Examples of systems with a discontinuous phase transition include 
underpotential deposition [121 Cll CH], while the adsorption of halides on 
Ag(lOO) [3 IHl 1201 mi are examples of systems with a continuous phase 
transition. 

4 Discontinuous Phase Transition 

A two-dimensional lattice gas with attractive adsorbate-adsorbate lateral in- 
teractions that cause a discontinuous phase transition is a simple model of 
electrochemical underpotential deposition [T21 CH UHl 1221 • Using a lattice- 
gas model with attractive interactions on an L x L lattice with L = 128, a 
family of FORCs were simulated, averaging over ten realizations for each re- 
versal curve at room temperature. The lateral interaction energy (restricted 
to nearest-neighbor) was taken to be (pij = 0nn = 55meV, where the posi- 
tive value indicates nearest-neighbor attraction. For this value of 0nn; room 
temperature corresponds to T = 0.8Tc, where is the critical temperature. 
The barriers for adsorption/desorption and diffusion (nearest-neighbor only) 
were Aa/d = Ann = 150 meV, corresponding to relatively slow diffusion P^ . 
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Simulation runs with faster diffusion (Ann = 125 meV) and the same adsorp- 
tion/desorption barrier showed httle difference from Fig. |21 indicating that 
diffusion effects are not significant for this model. The reversal electrochemi- 
cal potentials fir associated with the reversal curves were separated by 1 meV 
increments in the interval [— 200 meV, meV], and the field-sweep rate was 
constant at \dfi/dt\ = 0.3meV/MCSS. The FORCs are shown in Fig.H^a), 
with a vertical line indicating the position of the coexistence value of the 
electrochemical potential, fio = — llOmeV, and circles showing the position 
of the minimum of each FORC. 

In a simple Avrami's-law analysis, the FORC minima all lie at fl = fio 
However, in the simulations the minima are displaced. For 6 > 0.5, the 
minima occur at /i < /io, precisely at the points where the tendency to phase 
order, which drives local regions of the system toward the nearby metastable 
state {6 ~ 1), is momentarily balanced by the electrochemical potential, 
which drives the system toward the distant stable state ~ 0). For 6 < 0.5, 
the stable and metastable states are 6^1 and ^ 0, respectively, and the 
same balancing effect explains the FORC minima occurring at /i > /2o- 

The net effect is a 'back-bending' of the curve of minima, as seen in 
Fig. Efa). The definition in Eq. (Q) implies that the FORC distribution p 
should be negative in the vicinity of the back-bending. This can be seen 
in Fig|21^b), where the FORC distribution is plotted against the variables 
fib = {pT + /i)/2 and jlc = {fir ~ /^)/2- The negative values of p reflect a local 
divergence of the FORCs away from each other, which can be considered 
a dynamical instability, caused by the competition between the tendency to 
phase order and the effect of the electrochemical potential. It is important to 
note that, when the potential sweep is stopped suddenly, after the dynamical 
instability played out on short timescales, the system is observed to relax 
reliably to the stable state at that potential on large timescales. The only 
exception is the point {p = po, 6 = 0.5) along the FORC indicated by a 
bold line in Figs. Efa) andEl^b). It is also interesting to note that the curve 
connecting the minima of the FORCs resembles the van der Waals loop in the 
mean- field isotherm of a fluid system [^, but with an asymmetrical shape 
about the point {p = po, 9 = 0.5) and with a sweep-rate dependent shape 
(see Fig. El). 

We next explore this connection with the van der Waals loop in numerical 
solutions of a kinetic mean-field model. We can describe the competition 
between the phase-ordering and the influence of the potential explicitly, using 
a time-dependent mean-field model of the dynamics. The free energy is 
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Figure 2: (Color online.) First-order reversal curves (FORCs) for a discon- 
tinuous phase transition. The vertical line shows the position of the coex- 
istence value, /i = fLQ. The minima of each FORC are also shown (circles), 
(b) FORC diagram generated from the family of FORCs shown in (a). The 
positions of the FORC minima are also shown (circles). The straight line 
corresponds to the FORC for which the minimum lies at the coexistence 
value p, = JiQ (thick curve in (a)). After Ref. [T^ . 
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Figure 3: FORC minima dependence on sweep rate. The figure shows FORC 
minima for two families of FORCs with different sweep rates {fl = 0.03 and 
0.09 meV/MCSS). The lines are guides to the eye, obtained by smoothing the 
data using a Savitzky-Golay method [23 121] with a first-order polynomial 
and a window of 5 points. 



8 




Figure 4: (Color online.) Free energy of a fourth-order Ginzburg-Landau 
form as a function of 9, given by Eq. (0). The parameters were calculated as 
a = 30.2 meV, b = 16.0 meV, 9^ = 0.5, and {j2 - Jiq) = 0,2,4, and 6meV. A 
metastable state at 6* < 0.5 exists for the curves with (/i — /io) = 0,2,4 meV, 
but has disappeared in the curve {ft—flo) = 6 meV. The metastable minimum 
disappears at the so-called spinodal potential, /igp — /^o = Y\fi^ ~ 4.48 meV. 
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approximated by a fourth-order Ginzburg-Landau form |26j . 



F{e, m) = a ' - b'-^^ -{6- e^mt) - fio). (5) 

A plot of Eq. (jni) is shown in Fig. |3]for 6c = 0.5, with the values of a and b 
calculated as described below and with fi(t) — /xq = 0, 2, 4, and 6meV. In the 
noise-free (zero-temperature) case, the dynamics are given by 

^ = --^ = -- {a{e - Ocf - b{e - Oc) - m) - fio)) , (6) 

where 7 is a phenomenological damping parameter. The parameters to use in 
this mean-field model can be determined directly from the simulated family 
of FORCs. On the curve of minima of the FORCs, d6/dt = 0, which implies 

a{e - - b{e - e,) - - ^o) = o . (7) 



Differentiating with respect to /i and solving for dO/d/i gives 

d^ _ 1 

d^ ~ 3a{e - Oc)^ - b ■ 



(8) 



For 6 = 6c, this yields b = — (d^/d/x|6i=6»J ^- The so-called spinodal potential 
yUsp, and spinodal coverage ^sp, occur where d9/dfi diverges, so that 

3a(^sp - 0cf -6 = 0, (9) 

yielding 

a = 6/[3(^sp-^c)']. (10) 

The damping parameter 7 can be related to the "coercive potential" ficoer 
(The value of the electrochemical potential at which 6 = 6c on the major 
loop), the slope and the sweep rate dp,/dt as 



d^ 
d/2 



^ d9/dt\e=e, ^ (/icocr - /jp) 



where we have used Eq. (0) in the last step. 
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Figure 5: (Color online.) Reversal behavior of a dynamic mean-field 
model: (a) Families of FORCs. (b) FORC distribution for one family. 
The parameters were calculated as a = 30.2 meV, h = 16.0 meV, and 
7 = 3720MCSSmeV to fit the results of Fig. M^a). In (a), the sohd lines 
show the FORCs for faster sweep rate {Q = 0.3 meV/MCSS), and the dashed 
lines show the FORCs for a slower sweep rate {fl = 0.03 meV/MCSS). The 
thick line is the analytical result, Eq. ((Tj), for the line of minima, while the 
circles show the actual minima from the numerical integration. In (b), the 
full FORC distribution is plotted vs {fib, fie) (see text) for the FORCs with 
the faster sweep rate {Q = 0.3 meV/MCSS). The straight line corresponds 
to the FORC for which the minimum lies at the coexistence value ft = flo 
(thick curve in (a)). The appearance is similar to Fig|2fb). 
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Figure 6: (Color online.) KMC simulation results for FORCs with different 
sweep rates, = 0.03meV/MCSS (solid lines) and 0.09 meV/MCSS (dotted 
lines). The vertical line shows the coexistence value, p, = fto- The coercive 
potential for the FORCs with the slower sweep rate is shifted toward the 
fi = fio line. 



In Fig.El^a), we show the family of FORCs obtained by numerical integra- 
tion of the mean-field model, with parameters determined by a fit to the sim- 
ulated family of FORCs in Fig. Efa) . The line of minima in Fig. Efa) follows 
the analytical curve defined by Eq. (|7j), as expected. However, in contrast to 
Fig. 121^ a), the line of minima is symmetric about the point (/i = flo, 6 = 6c)- In 
addition, in KMC simulations at a slower sweep rate, we found that both the 
coercive potential (Fig. ^ and the line of minima (Fig. ^ of the FORCs were 
shifted toward p, = fio- In contrast, in the mean- field model at a slower sweep 
rate, the coercive field decreases, but the line of minima remains unchanged 
(as it must, since Eq. ((7j) is independent of the sweep rate). Preliminary sim- 
ulations indicate that including thermal noise in the mean-field model may 
further improve the agreement with the MC results. 
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Figure 7: (Color online.) (a) First-order reversal curves (FORCs) 
for a continuous phase transition simulated at a slow scan rate Q = 
0.0003 meV/MCSS. The thin black middle line shows the equilibrium curve. 
The inset is a magnification of the critical region. The minima of each FORC 
are also shown (black filled circles). The thick black line shows the first 
FORC which dips below the critical coverage, (b) FORC diagram generated 
from the FORCs shown in (a). The positions of the FORC minima are also 
shown (circles). The straight line corresponds to the first FORC for which 
the minimum dips below the critical coverage. After Ref. [T^ . 
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5 Continuous Phase Transition 



Using the same Hamiltonian, but with long-range repulsive interactions and 
nearest-neighbor exclusion as discussed in Ref. j7j, KMC simulations were 
used to produce the family of FORCs for a continuous phase transition. 
The reversal potentials ftr were separated by lOmeV increments in the in- 
terval [— 200 meV, 400 meV]. As in Ref. [2|, the repulsive interactions, 
with nearest-neighbor exclusion and (pnnn = — 21meV, are calculated with 
exact contributions for Vij < 3, and using a mean-field approximation for 
Tij > 3. The barriers for adsorption/desorption and nearest- and next- 
nearest-neighbor diffusion, are Aa/d = 300 meV, Ann = 100 meV, and Annn = 
200 meV, respectively [7j. Larger values of the diffusion barrier were also used 
to study the effect of diffusion on the dynamics. A continuous phase transi- 
tion occurs between a disordered state at low coverage and an ordered state 
at high coverage [201 1^ • The FORCs and the FORC diagram are shown in 
Fig. [7| Also indicated in Fig. U^sl) are the FORC minima and the equilibrium 
isotherm. 

Note that the FORC minima in Fig. [7|^a) lie directly on the equilibrium 
isotherm. This is because such a system has one stable state for any given 
value of the potential, as defined by the continuous equilibrium curve. The 
uniformly positive value of the FORC distribution in Fig. [7|^b) reflects the 
convergence of the family of FORCs with increasing /2. This convergence 
results from relaxation toward the equilibrium isotherm, at a rate which in- 
creases with the distance from equilibrium. It is interesting to note that, 
while it is difficult to see at this slow scan rate, the rate of approach to 
equilibrium decreases greatly along the first FORC that dips below the crit- 
ical coverage 6c ~ 0.36 (shown in bold in Fig. El^a)). The FORCs that lie 
completely in the range 6 > 0.36 never enter into the disordered phase, and 
thus their approach to equilibrium is not hindered by jamming. This is a 
phenomenon that occurs when further adsorption in a disordered adlayer is 
hindered by the nearest-neighbor exclusion. As a result, extra diffusion steps 
are needed to make room for the new adsorbates, and the system follows 
different dynamics than a system with an ordered adlayer |27]. The FORCs 
that dip below 6c = 0.36 enter into the disordered phase, and thus their ap- 
proach to equilibrium is delayed by jamming. This is reflected in the FORC 
diagram by the Florida-shaped "peninsula" centered around this FORC in 

Fig.ETb). 

The effect of jamming is more pronounced at higher scan rates, or with 
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a higher diffusion barrier, where the rate of adsorption is much faster than 
the rate of diffusion. The family of FORCs and FORC diagram at a higher 
scan rate, fl = 0.01 meV/MCSS, are shown in Fig. |H1 and the FORCs and 
FORC diagram with a larger diffusion barrier are shown in Fig. El In Fig. 
two distinct groups of FORCs undergoing jammed and unjammed dynamics 
can be clearly seen. This is reflected in the FORC diagram as a splitting 
of the "peninsula" into two "islands" of high p values. A similar effect is 
seen in Fig. |^ since also there the rate of adsorption is much faster than the 
rate of diffusion (larger diffusion barrier). However, Fig. El^a) shows a slight 
difference between the FORC minima and the equilibrium curve around the 
critical coverage. Notice also in Fig. |Hta) that even at a much higher scan 
rate than in Fig. [7| (nearly two orders of magnitude), the FORC minima still 
follow the equilibrium curve very accurately. Thus, the FORC method should 
be useful to obtain the equilibrium adsorption isotherm quite accurately in 
experimental systems with slow equilibration rates. 



6 Comparison and conclusions 

Two observations can be made by comparing the FORCs and FORC di- 
agrams for systems with discontinuous and continuous phase transitions. 
First, the FORC minima in systems with a continuous phase transition cor- 
respond to the equilibrium behavior, while they do not for systems with 
a discontinuous phase transition. Thus, FORCs can be used to recover the 
equilibrium behavior for systems with continuous phase transitions that need 
a long time to equilibrate. This could be useful in experiments. Second, due 
to the instability that exists in systems with a discontinuous phase transi- 
tion, the minima of the family of FORCs form a back-bending "van der Waals 
loop" , and the corresponding FORC diagram contains negative regions which 
do not exist for systems with a continuous phase transition. Since experimen- 
tal implementation of the FORC method should only require simple repro- 
gramming of a potentiostat designed to carry out a standard CV experiment, 
we believe the method can be of significant use in obtaining additional dy- 
namic as well as equilibrium information from such experiments for systems 
that exhibit electrochemical adsorption with related phase transitions. 
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Figure 8: (Color online.) (a) First-order reversal curves (FORCs) for a con- 
tinuous phase transition simulated at a high scan rate, = 0.01 meV/MCSS. 
The black curve in the middle shows the equilibrium isotherm. The minima 
of each FORC are also shown (black filled circles), (b) FORC distribution 
generated from the FORCs shown in (a). The positions of the FORC minima 
are also shown (black filled circles). The straight line corresponds to the 
FORC for which the minimum lies closest to the critical coverage. 
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Figure 9: (Color online.) (a) First-order reversal curves (FORCs) for a con- 
tinuous phase transition simulated with a large diffusion barrier Ann = 300 
meV. The black curve in the middle shows the equilibrium isotherm. The 
minima of each FORC are also shown (black filled circles), (b) FORC distri- 
bution generated from the FORCs shown in (a). The positions of the FORC 
minima are also shown (black filled circles). The straight line corresponds to 
the FORC for which the minimum lies closest to the critical coverage. 
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Appendix 

In this appendix we present a mapping between lattice-gas models of adsorp- 
tion and discrete spin models of magnetic systems, and then introduce the 
FORC method in the original magnetic language. 

The occupation variable in the lattice-gas model, q G {0, 1} is a binary 
variable, just like the magnetization variables: Si = Ms/rris G {—1,1} in 
the Classical Preisach Model (CPM). We therefore have the mappings Cj = 
(si + l)/2 and {fi- fio) = 2H As a result, the FORC method can 

be apphed to electrochemical adsorption, as well as to magnetic hysteresis. 

The CPM is based on the idea that a material consists of a number of 
elementary interacting "particles" or "domains," called hysterons. The hys- 
terons are assumed to have rectangular hysteresis loops between two states 
that have the same magnetization values, +ms and — m^, for all hysterons. 
A typical hysteresis loop for a hysteron is shown in Fig. ^| Hu and Hd are 
the up and down switching magnetic fields respectively. It is also assumed 
that the different hysterons have a distribution of reversal fields ^{Hu, Hd). 
In the CPM, the total magnetization can be defined as [30] 

M{t)= ff <l>{H^,Hd)[m{Hu,Hd)H{t)]dH^dHd , (12) 

where the operator m{Hu, Hd) applied to H{t) gives +ms if the particle is 
switched up and — if the particle is switched down. Note that rh{Hu, Hd)H{t) 
depends on the field history of H{t) and not only on the instantaneous value, 
which enables the CPM to model irreversible hysteresis behavior. For a de- 
tailed discussion of the Preisach Model, see Ref. pUj . 

In a typical FORC analysis of a magnetic system, the magnetization M is 
saturated in a positive applied magnetic field, and then the applied magnetic 
field is decreased continuously to a reversal field Hj,. The magnetic field is 
then increased back to saturation. A first-order reversal curve is the response 
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Figure 10: Schematic diagram of a hysteresis loop for a single hysteron. if, 
and Hd are the up and down reversal fields, respectively. After Ref. [HUj . 



of the magnetization to the increasing magnetic field {H > Hr). This is done 
for different values of Hr, and a set of curves, M{Hr,H), is collected PI ITT]. 
The FORC distribution is defined as 

. 1 d'MjHr, H) 

where the tilde denotes a trivially different normalization from the one used 
here. Thus, using the mapping given above, one arrives at Eq. as the 
definition of the FORC distribution in an electrochemical system. 
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